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1/-^ ■ Abstract 

Q ' The thermodynamics of vortices in the critically coupled abelian Higgs model, defined on 
^ ■ the plane, are investigated by placing N vortices in a region of the plane with periodic boundary 
^ • conditions: a torus. It is noted that the moduli space for vortices, which is the same as that 
of N indistinguishable points on a torus, fibrates into a (DP/v_i bundle over the Jacobi manifold 
^ of the torus. The volume of the moduli space is a product of the area of the base of this bundle 
and the volume of the fibre. These two values are determined by considering two 2-surfaces 
J> ' in the bundle corresponding to a rigid motion of a vortex configuration, and a motion around 
^ a fixed centre of mass. The partition function for the vortices is proportional to the volume 
■ of the moduli space, and the equation of state for the vortices is P{A — AnN) = NT in the 
thermodynamic limit, where P is the pressure, A the area of the region of the plane occupied 
by the vortices, and T the temperature. There is no phase transition. 
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1 Introduction 



Solitons, as topological defects in field theories and models for elementary particles and mag- 
netic flux tubes in superconductors, have been attracting a lot of interest recently. Defects 
such as strings and monopoles, if they exist, are recognised to have a significant effect on the 
dynamics of the early universe 0]. Solitons in Skyrme's model, a low energy limit of QCD, are 
a good model for baryons and their low energy interactions 0. 

In general, dealing with the interactions of solitons requires a treatment of the full field 
theory in a non-perturbative regime. However, in certain theories at low energies, the motion of 
a soliton configuration can be approximated by a continuous sequence of static configurations 
0, in much the same way that a moving cinema picture is built from a sequence of still 
pictures. The energy functional of the fields stays "close" to its minimum, which occurs on the 
configuration space of stable static solutions of the field equations. This approximation rules 
out higher energy phenomena, like soliton-antisoliton pair creation, as there are no stable static 
configurations containing a mixture of solitons and antisolitons. The finite dimensional manifold 
parametrising static configurations is referred to as the moduli space, and the interaction of 
the solitons in the field theory is modelled by geodesic motion in this space, with a metric 
determined from the field theory Lagrangian. This approximation is expected to be good for 
low soliton impact velocities, becoming exact in the limit of soliton velocities tending to zero. 

The moduli space approach to sohton dynamics is possible at a critical value of a coupling 
constant in the Lagrangian, where the equations describing the minima of the energy functional 
reduce to first order differential equations, the Bogomolny equations |Q, and there are no forces 
between the solitons. This is why the moduli space has no potential, and the motion on it 
is geodesic motion. Comparison of the geodesic motion with numerical simulations of the 
scattering of critically coupled vortices shows a good agreement for velocities up to around |c. 

For many applications, it would be useful to know the thermodynamic behaviour of a 
large number of solitons. In a recent paper by one of us [^], it was shown how to derive the 
thermodynamics of vortices in the abelian Higgs model at critical coupling, with the vortices 
moving on a 2-sphere. The calculations were performed analytically, but within the framework 
of the moduli space approximation. This method can be applied to vortices moving on a general 
orientable compact Riemann surface S. The Lagrangian density of this model is 

- + lD,<p{Dyr - Im' - ir a) 

(j) is complex scalar Higgs field and D^cp = d^cj) — ia^cp ; where is the U{1) gauge potential. 
ffiu = d^ay—dutt^ is the Maxwell field tensor. Although the Bogomolny equations describing the 
fields of minimal potential energy are first order, explicit solutions are not available. However, 
a static solution with energy A^vr and total magnetic flux 

Jh2 = 2nN , Nez (2) 

is uniquely determined (up to a gauge transformation) by specifying the positions of the pre- 
cisely zeros of the Higgs field, which may be identified with the positions of A^ vortices 
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Classically, these vortex positions are an unlabelled set; the configuration is unaffected by inter- 
changing any of the vortices so the moduli space A^at is (S)^ /T,]\f where Sat is the permutation 
group of N elements. 

There is a naturally defined metric on A4n, derived from the kinetic part of the Lagrangian. 
Let QijiQ) denote this metric, where Q = {Q^} is a set of arbitrary coordinates. In the moduli 
space approximation the energy of moving vortices is \T^gij{Q)Q^Q^ ■ By introducing conjugate 
momenta Pi = T^gij{Q)Q-' ■, the energy may be written in the form 

E{P,Q) = ^g^^{Q)P,P, (3) 

The statistical mechanics of vortices at temperature T can be treated by using a Gibbs 
distribution, where the partition function is 

Z = J^JdP){dQ)e^^ (4) 

and h is Planck's constant. By doing the Gaussian momentum integrals, this reduces to 

Z={f {dQ)^d^^{^r (5) 
Jmn ^ h 

The first factor is simply the volume of the moduli space. It is known that the metric on the 
moduli space is Kahler and hence the volume of the space can be determined from the areas of 
2-surfaces in it by using homology and symmetry arguments. 

In this paper, we investigate the statistical mechanics of vortices on a torus rather than on a 
sphere, thereby eliminating undesirable curvature effects. The torus corresponds to a unit cell 
in the plane with periodic boundary conditions. One expects that in the thermodynamic 
limit — > oo at fixed number density n = ^, the behaviour of the particles should not depend 
on the global topology of the manifold they are defined on, although this has not been proven. 



2 The fibre bundle structure of the moduh space 

To model a torus M = 5*^ x 5^, we choose a unit rectangular cell in the plane with periodic 
coordinates (xi, X2) with ranges Xi G [0, 1), X2 G [0, a), and give the spacetime the metric 

ds^ = dxl - L'^idxf + dxl) (6) 

where L is a constant physical scaling factor; the area of the torus is then A = L'^a. We 
introduce a complex coordinate z = xi + 1x2 ■ The Lagrangian (|ID leads to the Bogomolny 
equations for critically coupled vortices 

{Di + iD2)(j) = (7) 

fi2 + Il\\<I>\' - 1) = 
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Here, and are not strictly speaking functions on the torus, but rather a section and 
connection on a U{1) bundle over it. The total flux is quantised as in (H), and is the first 
Chern number of the bundle. 

The moduli space for N vortices on the torus M is Mn = M^/Stv- This makes it appro- 
priate to define the configuration in terms of its divisor d{(f)). A divisor of a function is a formal 
unordered sum of the zeros (counted positively by multiplicity) and poles (counted negatively) 
of the function. The divisor of the Higgs field of N vortices at positions {pi, . . . ,Pn} is 

N 

di^) = Y.l-Pi (8) 

i=l 

The magnitude of the Higgs field away from the vortex positions is not important for our 
purposes. Let us therefore consider the complexification of the U{1) gauge group, namely (D*. 
Under transformations of the form e*°'-^'^V 5 where a{z, z) G (D, the absolute magnitude |0| 
is no longer gauge invariant but the positions of the zeros of 0, and hence the vortex positions, 
still are. Using the enlarged gauge group, we can locally put the fields in the holomorphic 
gauge where = |(ai + ia2) = 0. Then the first Bogomolny equation, which can be written 
as {dz — iaz)(p = 0, reduces to 

I- 

so becomes an analytic function of z. In this gauge we can express the Higgs field as a product 
of analytic functions with precisely one zero each on the torus. These functions are the theta 
functions, which are actually defined on the universal covering space of the torus, the complex 
plane (D. The torus can be represented as the quotient of (D by the lattice generated by 1 and 
T = ia. The theta functions are not singled valued on this quotient, but the positions of the 
zeros in the complex plane map to a single point when quotiented by the lattice. 

The lattice generated by 1 and r defines a normalised choice of what is known as the Jacobi 
mapping. The Jacobi mapping allows us to define a centre of mass for the vortex configuration 
(which is not unambiguously possible on the original torus). It is this property that allows a 
fibration of the moduli space into coordinates related to (but not the same as) the centre of 
mass and relative positions of the vortices. This uses the theory of generalised theta functions, 
quasiperiodic functions (or complex analytic relatively automorphic functions" - carafs for 
short) which can be defined for any orientable Riemann surface of genus g > 0. Our discussion 
of them will be brief and specific to the torus, but the interested reader is referred to [0] for a 
fuller and more general discussion. 

The Jacobi mapping of an orientable Riemann surface S of genus g maps 5 to a complex 
^f-torus, and is defined in terms of the integrals of normalised abelian differentials (of the first 
kind). On the torus, there is just one normalised abelian differential, namely dz, whose integral 
is z, so the Jacobi mapping is essentially the identity. The fundamental periods of dz are 1 
and T. More precisely, let us choose a base point po on M, and identify it with the origin of 
coordinates. We define the Jacobi manifold of the torus J(M) by 

C 

J(M) = - , £ = {ni + nsT : rii, G Z} (10) 
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and the Jacobi mapping $* by 



$* : M J(M) 

^*{xi,X2) = Xi+ iX2 



xi e [0, 1),X2 e [0,a) 

Given the Jacobi mapping there is a natural map from an A^-vortex configuration to (D defined 
by = I^ili ^*iPi)- ji is then a particular choice of a centre of mass for the configuration. 
The simplest example of a 'classical' theta function is 

n=+oo 



-i (-l)"e("+2 

n=— oo 



(12) 



which has zeros at the origin in the complex plane and lattice translates of that. Translation of 
this function by lattice vectors has the effect of multiplying it by non- vanishing factors, termed 
factors of automorphy: 



Oi{t,z + 1) 



-Oi{t,z) 
-k{t, z)9i{t,z) 



(13) 



-27Tiz ^ — TTIT 



'we shall now drop the r argument and just write 6i{z)). There are 



where ^(r, z) = 
three other classical theta functions, related to 6i {z) by 



9^{z) 



-d2iz + -) = -ie 



-e''''esiz + - + -T) 



-ie—e'^^'e^iz + -t) 



(14) 



Various identities exist between these four functions; the principal ones that we will use are 
summarised below: 





e,{z) 


02{Z) 




0,iz) 


1 


-1 


-1 


1 


1 


r 


— K 


K, 


K 


— K 



(15) 



The factors of automorphy of classical theta functions 



om = 02{-) 

el{z)el 
el{z)el 

Oiiy + z)ei{y - z)el 



el{z)el - el{z)el 
6l{z)el - el{z)el 
el{y)el{z) - el{y)el{z) 



(16) 



where 6i denotes 6i{0). The reader is referred to, for example for an extensive list of the 
properties these functions possess. These functions can be used to represent divisors on the 
torus M. Although they are not single- valued on M, the factor of automorphy is non- vanishing, 
and hence the zeros of 6i{z) in the complex plane map to a single point of M. 
In the holomorphic gauge we can write the Higgs field locally as 



N 



(17) 



i=l 
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which is an example of a caraf. Under translation by a lattice element A, this has factor of 
automorphy 



(-1)^ for A = 1 

UiLi z - Zi) for A = r 



The characteristic c(^) of a factor of automorphy ^ of a scalar caraf is equal to the order of the 
divisor of the function, i.e. the number of zeros minus the number of poles. Hence 0(^0) = A^. 
It is a special case of a general theorem on factors of automorphy that we can write 



e^(A,z) = p,(A)C(A,^) (19) 

where t = t^cj)) = Y^iLi Zi as defined above, denotes the factor of automorphy of Oi{z)^ and 

pt{l) = l pt(r) =exp27r2t (20) 

By another theorem, the carafs on the torus with factor of automorphy Pt^C^ form at fixed t 
a vector space with dimension A^. Any basis functions of this space, /j(t, z) {i = 1 . . . N) are 
quasiperiodic in t as well as z: 

Mt,z + X) = p,{X)^^^{X,z)Mt,z) (21) 
Mt + X,z) = p,{X)x^J{X,t)f,{t,z) (22) 

where p2(A) is defined identically to pt(A). The vector valued function 0{t,z) = {fi(t,z)) is 
termed a generalised theta function of rank A^, with associated theta factor Xiji^^ t)- 

We now give a useful basis for these functions on the torus for general rank A^, and give the 
theta factors for first and second rank. For rank 1, a basis having the factors of automorphy 
(^ and (0) is 

fit,z) = e,iz-t) (23) 
Using the properties of Oi{z) listed in (^3|) the theta factor can be computed to be 



X(l,t) = -1 (24) 
x{T,t) = -e-'^^e-^^'' = -K{T,t) 

For rank 2, a basis is provided by the independent functions 

f,{t,z) = {el{z-^-),el{z-^-)} (25) 

The computation of the theta factor utilises identities on the squares of theta functions, and 
leads to 

Xil,t) = %] (26) 



X{r,t) 



-k{t, I) 
-ft:(r, I) 



' 2' 
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For rank even, we will find the following basis useful 



f,{t,z) = {e^,e^-'el...,e^, (27) 

where all these theta functions are evaluated at 2; — For odd we shall use 

f,{t,z) = {9^,9^-'el...,e,e^-\ (28) 

Ol (74^2^3, 04^^203, ■ ■ ■ , [71^4 (72(73, 174 (72^3/ 

These functions all have the same factors of automorphy, and are clearly independent, as they 
have zeros of orders A^, A^ — 2,A^ — 3,...,2,1,0 at z = and lattice translates of it. It is 
interesting to note the absence of a basis element with a zero of order A^ — 1. If such an 
independent function existed, division by 0^{z — would lead to an elliptic function with one 
simple pole in a unit cell, but this contradicts a theorem on elliptic functions. Computation of 
the theta factors of these bases is in principle possible, but in practice difficult; we shall anyway 
not need to use the theta factors in our computation of the volume of the moduli space. 

The relevance of these generalised theta functions and their bases to the moduli space of A^ 
vortices is clear from the following theorem |^ : 

Theorem 1 For an orientable Riemann surface S of genus g, and any integer N > 2g ~ 1, 
the space — — {d^ : d^ — X]i=i ^ ' Pi ; Pi ^ of positive divisors of order N on the surface 
can be given the structure of a fibre bundle with fibre (DPat-^ and base space J{S), the Jacobi 
manifold of S . This bundle has a natural projection 

$^ : ^-J(5) (29) 
^%ipi...pN) = $*(pi) + ■■■ + $* (pjv) mod £ 

where C is the lattice of periods of normalised abelian differentials. 

We can understand this for the torus M, with g = 1, a.s follows. $|^r may be identified with 
t mod £, with t defined as above. As noted above, this coordinate is related to a choice of 
centre of mass of a configuration; its most important property is that it is invariant under the 
permutation group acting on the vortices. For fixed t we have the A^-dimensional vector space 
with a basis given by (^) or (|28|) above, and a divisor determines an element of this vector 
space via the intermediary of the function 

N 

Mu, t,z) =Y^ Uifi{t, z) (30) 

i=l 

such that d{f^{u,t, z)) = d^f, for a choice of m,. The Ui give homogeneous coordinates in the 
(DPat-i fibre, and are only defined up to a constant (non-zero) complex multiple, as cui rep- 
resents the same divisor as Uj. If t is restricted to lie in one unit cell of the lattice, the 
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Ui are unique modulo such a constant multiple. The t coordinate is unique modulo lattice 
translations, and translation of t by a lattice spacing changes the coordinates Ui by the appro- 
priate thcta factor (dependent on the choice of caraf basis functions fi{t,z)). More explicitly, 
(t, Ui) (t + A, UjXji^{t, A)) (where an irrelevant homogeneous factor of p^^{X) has been omit- 
ted) leads to the same divisor, as writing t' — t + X and u'^ — Ujxjl 

U{u',t',z) = ^w^/,(t + A,^) (31) 

= PzWU{u,t,z) 

Hence d{f^{u', t' , z)) — d{f^{u, t, z)). The theta factor therefore describes the non-trivial nature 
of the bundle over J(M). 

With the understanding that t will now be restricted to lie in the unit cell containing the 
origin in the complex plane, we now seek to investigate the form of the metric on the bundle, 
using the symmetry of the torus. 

3 The metric on the moduh space M.n 

We begin with the most general Hermitian metric on the moduh space, which has the structure 
of a fibre bundle as discussed in section 2 : 

ds^ — a{t, v)dtdi + ba{t, v)dtdVa + ^a(^, v)dtdVa + Ca/3(i, v)dvadvj3 (32) 

Va are the inhomogeneous coordinates on the fibre CPat-i, a, f3 — 1 . . . N — 1 , and t is the 
base coordinate which covers J{M) once, a is real and Ca/s is Hermitian. We note first that 
from translation symmetry of the torus, a,ba, Cap must be functions of v only. Furthermore, we 
can use the 180° rotational symmetry of the torus to eliminate the functions ba- The proper 
distance associated with a small displacement {St, Sv) satisfies 

Ss'^it, V, 6t, 6v) = Ss\-t, V, -6t, 6v) (33) 

and by choosing St real so that St — St, this implies that 

bJVa + baSVa = (34) 

Choosing Sva — iSap implies that Im{bp) — 0, and choosing Sva — Sap implies that Re{bp) — 0, 
hence 6/3 = for each /3. 

It would seem physically reasonable that for a small rigid motion of the vortex configuration, 
which would correspond to shifting the t coordinate but staying stationary in the fibre, the 
distance moved in the moduli space should not depend on the relative positions of the vortices. 
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This is confirmed when we apply the resuh that the moduh space is Kahler. The Kahler form 
is 

uo = -{a{v)dt A dt + Cai3{v)dva A dvfs) (35) 

On a Kahler manifold du = , which implies 

da J 7 7- da , , dcas , , dcaa . , ^ /o^\ 

— — dVa A at A at + -r^dVa Adt Adt + — — dv^ A dv^ A dvg + ^_ dv^ A dVa A dvn = (36) 

OVa OVa OV^ OVj 

SO 

and therefore a is a constant on the fibre bundle. The resultant form of the metric is then 

ds'^ = adtdt + Caf3{v)dvadvf3 (38) 

so the volume of the moduli space is a product of the area of the base J{M) and the volume 
of the fibre (DP^_i. 

Samols has given a formula for the metric at a point on the moduli space in terms of the 
corresponding solution of the field equations. Writing / = log |0p (in the original formulation 
with gauge group U{1)), it follows from the Bogolmolny equations (0) that / satisfies 

N 

Q2f + L\1- ef) = 47r ^ 5^^\z - y,) (39) 

i=l 

where Ui are the vortex positions in J{M). f can be expanded as a series about each vortex 
position 

f{z,z) = \og\z-yi\'^ + ai{y) + ^bi{y){z-yi) + h)i{y){z-yi)+ (40) 
Ci{y){z - yif + Ci{y){z - y^f + di{y){z - yi){z -?/,) + ... 

where y denotes the set of all vortex positions, bi measures how the centres of the contours of 
|0p drift away from the vortex position yi as the magnitude of |0p increases from zero, and 
this depends on the positions of all the other vortices. The metric on the moduli space is then 

^ Bi- 
ds' = Y: {L% + 2-^)dydm (41) 

i,j=l (^Vi 

and we shall show in the next section how to relate this to 



4 Collective motion of N coincident vortices 

We now turn to the calculation of the first factor in the volume oi M.n '■ the area of the base of 
the fibre bundle. From (|5BD, this is aa, since the range of t is {0 < Re{t) < 1 , < Im{t) < a}. 
To calculate a we look for a 2-surface of vortex configurations that is locally orthogonal to the 
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fibre. Pursuing our observation that the base coordinate defines a choice for the centre of mass 
of the vortices, such a surface should correspond to a rigid motion of the vortices. A simple 
choice of such a surface is obtained by placing all the N vortices at one point, and allowing 
that point to vary around the torus. We define 

Mc = : = N-p,peM}c-^ (42) 

We can express Aic algebraically in terms of the fibre bundle coordinates as follows. Write the 
Higgs field in the holomorphic gauge as 

0(^) = <(^-y) (43) 

where y = $*(p) G J{M). Recall the basis of carafs of rank given in (^) and (pSf) , 
fi = {9i{z — j^), . . .}, where now t = Ny mod C. Keeping y restricted to the small patch of 
the unit cell delimited by < Re{y) < , < Im{y) < , which we shall denote lAn (we 
shall use Urs-, r,s = 1 ... to mean the small patch ^ < Re{y) < j^, -^^^^ < Ira{y) < ^ ), 
then 

0(^) = /i(Ary,z) (44) 
Hence the coordinates of A^c for y E lAn are 

Ma = {(t, u) : {Ny, 1, . . . 0) , t G J(M) , u E CPjv-i} (45) 

This surface is locally orthogonal to the fibre since u is constant, and as y varies over Uu, t 
varies over the whole of J{M). To find the coordinates of the surface for y E Urs one must 
apply the quasiperiodicity of the caraf basis in the t coordinate, corresponding to the non-trivial 
structure of the fibre bundle. For y varying over the whole of J{M), the surface intersects each 
fibre A^^ times and always orthogonally. For example, for y E Uu, equation (p2|) implies 

u^iy + l^) = Xij{hNy)uj{y) (46) 

where the LHS are the coordinates in U21, and a homogenous factor has been omitted. Using 
A^ = 2 as an example, we have determined the xi^^t) matrices and hence the surface Aic has 
coordinates 

(2y,l,0) yEUu 
(2y- 1,^1,^2) yE U21 
(2y-r,0,l) yEUu 
{2y-l-r,elel) yE U22 

Returning to A^ vortices, the algebraic expression for the coordinates of A^c can be inserted in 
the metric on the fibre bundle, restricted to A4c- For y E Uu, dt = Ndy and dv = ^dy = so 
recalling equation ( ^8l) 

ds^ \m,= aN'^dydy (47) 
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We can now find a by determining the metric on Ai^ from Samols' formula (^) which for 
coincident vortices reduces to 

db 

ds^ \m.= NiL^ + 2—)dydy (48) 

oy 

where ^b{y) is the hnear coefficient in the expansion of / = log |0p around the vortex position 
y. The reflection symmetry z — y ^ —{z — y) can be used in this local expansion, and implies 
b{y) = ; the contours of the Higgs field are locally centred on the A^-vortex position. Hence 
(^81) simplifies to 

ds'^ \j^^= NL^dydy (49) 

which is the same result as would be obtained for coincident vortices in a plane, using the 
more powerful circular symmetry. Comparison with (|47|) yields 



N 



(50) 



Hence 

Area{J{M)) = ^ = ^ (51) 

5 Simultaneous motion of two diametrically opposite 
vortices 

In the previous section, we used a rigid collective motion of the vortex configuration to determine 
the area of the base space of the fibre bundle. We now turn to determining the volume of the 
fibre, (DP/v_i, by looking for a motion that lies entirely in the fibre, keeping the centre of mass 
fixed. The simplest example seems to be the motion where — 2 of the vortices are fixed at the 
point po, and the other two are at positions pi and p2 symmetrically placed on opposite sides 
of po and moving over some suitable region of the torus. This motion of vortices corresponds 
to a motion in a submanifold M.d of the moduli space, 

Md = {d{(t)) : rf(0) = (AT - 2) • po + 1 ■ Pi + 1 ■ P2, Pi, P2 e M} (52) 

where $*(pi) = x, say, and $*(p2) = mod C , so the coordinate t(0) = 0. M.d is covered 
once if x ranges over the half cell Wi = {s : < Re{x) < |,0 < Im{x) < a}. As before, we 
seek to find an expression for the fibre coordinates on A^^ . Again in the holomorphic gauge 

(j){z) = ei{z - x)ei{z + x)e^-\z) (53) 



We can use the addition identities of theta functions ([T6|) to write this as 

0W = ^«(^)^4(^) - el{z)e^''{z)el{x)) (54) 
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Recall our basis, with t — here, 



Mo,z)^{e^{z),e^-\z)el{z),...} (55) 

Hence, omitting the irrelevant 9^^ factor, the base and homogenous fibre coordinates of Aid 
are 

Md = {{t, Ui) : (0, el{x), -el(x), ... 0)} (56) 

Therefore Ma lies in the CPi line U3 — U4 — . . . — un — 0. The inhomogenous fibre coordinates 
are 

v^=^ eJM " - ^ (57) 
otherwise 

which are even elliptic functions of x, as expected from the symmetry of the motion. Now for 
X varying over a unit cell, the number of solutions of i;i = c G C is independent of c, and is 
equal to the number of zeros or poles of Vi, counted by multiplicity. In this case, there are two 
solutions, but the symmetry Vi{x) = Vi{—x) implies that there is exactly one solution in the 
half cell Ui. Hence, for x varying over Ui, Aid covers a (DPi line of (DPat-i exactly once. 

We now need to apply this result to the metric on the fibre bundle, whose restriction to Add 
is given by 

ds^lj^ =2(L'^ + ^ -^^)dxdx (58) 
ox ox 

(remembering that as we vary x we are varying the positions of two vortices simultaneously). 
^bx is the linear coefficient of the expansion of / = log|0^| around the vortex at x and ^b-x 
the linear coefficient around the vortex at —x. Set bx = b, and by 180° rotational symmetry 
b-x = —b. However, it is now not possible to solve for b using the symmetry of the torus. 
Fortunately, we do not need to know b everywhere; we can instead consider the integral of 
over X eUi, and use Stokes' theorem. The area of Add is 
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7 = / (2L^ + 2^-2^){-dxAdx) (59) 
iw/ dx dx '^2i ' ^ ' 

\% / dx{bdx) 

JlAi 



2 

= A-2i 

k 

s 

= A + 2i I bdx 

JdUi 

2 



where the orientation of the contour dUi is shown in figure 1. b has poles at the points where 
the vortices at x and — x coincide, and these are marked. Consider the integral on the two line 
segments h = {x : x = ct,c e (e, | — e) U (| + e, 1 — e)}. For x e h, \(j>{z)\ = \(l){—z)\ = \(j>{z)\ = 
|0(-^)|, and 

' bdx = (60) 
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as the parts from c < | and c > \ cancel. In a similar way, the integral on 1^ (see figure) is 
zero, and I2 and Z4 are lattice translations of each other with opposite direction, and so their 
contributions to / cancel. The only contribution to / comes from the poles indicated. Consider 
the pole at x = 0. To investigate it, we expand / about ^ = 0, assuming x is near 0. 

/ = {N -2)\og\z\'^ + \og\z + x\^ + \og\z-x\'^ + A{x)+ (61) 
-B{x)z + -B{x)z + C{x)z^ + C{x)z'^ + D{x)zz + . . . 

where A,C,D are all regular around x — The symmetry z ^ —z implies that B{x), and 
all other coefficients of terms odd in 2;, vanish. Likewise, the symmetry x —x implies 
A{x), C{x), D{x) are all even functions of x. To find the singularity in 6, compare this with an 
expansion about z = x 

1 1- 

/ = \og\z - x\'^ + a{x) + -h{x) {z - x) + -b{x) {z - x) + (62) 
c{x) {z — xY + c{x) {z — xY + d{x) {z — x){z — x) + . . . 

Using log \z + = log \2x\^ + ^{z - x) + ^{z - x) - ^{z - xf - ^{z - xf + . . . (the 
expansion is valid for 1^; — a;| < \2x\, which includes z — 0), and similarily for log we obtain 

a{x) = A{x)+\og\2x\'^ + {N-2)log\x\'^ + C{x)x'^ + C{x)x^ + D{x)xx + 0{x^) (63) 
^ b{x) = — + — — - + C{x)2x + D{x)x + 0{x^) 



2 ^ ' 2x X 

1 N -2 

8a;^ 2x 
d{x) = D{x) + 0{x'^) 



Hence, a{x) , b{x) , c{x) have logarithmic, first order and second order poles at x = 0, with the 
singularities completely determined by the log terms in the expansion around z = 0. The 
integral of 2ib on the contour near x = 0, = {x : x = ee*^, 6 G (0, ^)} just has a contribution 
from the pole (the linear and all higher order terms are odd in x and vanish under integration): 



2i bdx = 2i (2N - 3)ide (64) 

Jle Jo 

= -7r(2iV-3) 

A calculation for the pole of 6 at x = |, where the two vortices again coincide (but this time the 
expansions are unaffected by the other vortices at z = 0), yields a contribution to the integral 
of — 27r (this time the integral is around a half circle), and likewise for the pole at a; = \ + \- 
The poles at a; = |, | + r contribute — tt as the integral there is around a quarter circle. Finally, 
the pole at a; = r makes the same contribution as a; = 0. The residues of all the poles of h on 
dlAi are shown in figure 1. Their contributions to the integral / sum to give 

I = A-27r(2iV-3 + 2 + l) (65) 
= A-AtiN 
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The result is always positive, because of the observation in that the Bogomolny equations 
for vortices on a compact Riemann surface only have non-trivial solutions if the physical area 
A of the surface is larger than AtcN. 

The value of / obtained for x eUi represents the area of a complex line, (DPi, in the fibre 
(DPiv-i (had we integrated instead over the whole unit cell U, we would have obtained double 
the result from the sum of pole residues, but would then divide by two as this surface would 
have covered CPi twice). Now in a complex manifold with a Kahler metric, the area of any 
complex curve is the integral of the Kahler form over it. So the integral of the Kahler form 
(|35D over a complex line in the fibre is 



U = / uj = A- AnN (66) 

JCPi 

u = u!/[uj] has integral of unity over a complex line, and so is a basis element for H2{(CPn-i, Z), 
the second cohomology group of the fibre over the integers. Then, from the topology of complex 
projective space (see e.g. |1^), u}^~'^ is a basis element for H2n-2{^Pn-i,'^)- Now the volume 
of a Kahler manifold of dimension 2N — 2 is the integral of uj'^~^/{N — 1)! over it. Therefore 
the fibre volume is 

6 The partition function of N vortices 

We can now easily combine our results for the area of the base space of the fibre bundle (|5T|) 



and the volume of the (DP/v_i fibre (|67|); their product gives the volume of the moduli space 



A 

VoI{Mm) = -j^{A - 47riV)^-i = ^(1 - 4™)^-i (68) 
where n = ^ is the vortex number density. The partition function is then 

Ztorus = ^{l-^7.nr-\^-^r (69) 
This should be compared with the result obtained in 

Zs,Here = ^il-^7rnrC-y-r (70) 

The additional single factor of (1 — Aim) is thermodynamically insignificant. We can calculate 
the free energy F = —T In Ztoms, using A^! ~ In — for large A^, 

2e7r^T 

F ~ -A^T(-lnA^ + ln(A-47rA^) + ln— — ) (71) 



The pressure P = — 14- is 



dA 

NT 

P = —-^ (72) 

A - AnN ^ ' 
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and the entropy S 



is proportional to A^: 



S = N{- In + In (1 - 4?™) + In 



2eW 



(73) 



So, in the thermodynamic hmit — *^ oo at fixed n, we obtain the same physical results on the 
torus as on the sphere, a remarkable result given the disparity of the mathematical methods 
used to arrive at the result for the two manifolds. This result supports the general principle that 
in the thermodynamic limit the behaviour of the gas is independent of the global topology of 
the physical manifold it occupies, even though topological methods have been used to calculate 
the partition function. Since A — > oo as — oo and since the torus is intrinsically flat, we can 
interpret the results ([72D and (^) as describing the thermodynamics of vortices in the plane. 
In the low density limit we flnd 



as would be obtained for a gas of hard discs of area 27r, and in the high density limit n = 
with e small 



This describes what would happen if, keeping the number of vortices flxed, we were to 
decrease the physical area available to the vortices towards the limit A = 4ttN; the work 
needed to compress the surface — / PdA would diverge as A — AnN. 

We also note that the pressure and free energy are differentiable, provided A > AnN, and 
hence there is no phase transition in a vortex gas at critical coupling. This is perhaps not 
unexpected due to the smooth interactions of the vortices and the absence of forces between 
them. 

It would be interesting to investigate the behaviour of the vortices some small distance away 
from the critical coupling either towards a Type I superconductor (so that the vortices would 
be weakly attracting) or a Type II superconductor (with vortices repelling). In these cases, 
it is possible that the vortex gas could undergo a phase transition between a gaseous and a 
condensed (type I) or crystalline (type II) phase. 



PA = NT{1 + Ann + 0(n^)) 



(74) 



P = Te-i + 0(1) 



(75) 



15 



References 

[1] E. Kolb and M. Turner, The Early Universe (Chicago University Press 1990) 

[2] T.H.R. Skyrme, Nucl. Phys. 31(1962) 556 ; E.M.Nyman and D.O.Riska, Rep. Prog. Phys. 
9(1990) 1137 

[3] N.S. Manton, Phys. Lett. B110(1982) 54 

[4] E.B. Bogomolny, Sov. J. Nucl. Phys. 24(1976) 449 

[5] N.S. Manton, Nucl. Phys. B400(1993) 624 

[6] S. Bradlow, Commun. Math. Phys. 135(1990) 1 

[7] R.Gunning, Riemann Surfaces and Generalised Theta Functions (Springer Verlag 1976) 

[8] E.T.Whittaker and G.N.Watson, A Course of Modern Analysis (Cambridge University 
Press 1980) 

[9] T.M. Samols, Commun. Math. Phys. 145(1992) 149 ; PhD Thesis, Cambridge University 
1990 (unpubhshed) 

[10] P. Griffiths and J. Harris, Principles of Algebraic Geometry (Wiley 1978) 



16 



_ 2N-3 1. 1 




2N-3 h 1 x=l 

Figure 1: Integration contour dUi and residues of b 
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